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ABSTRACT 

We present Spitzer IRAC imaging of the large-scale jet in the quasar PKS 1136-135 at wavelengths of 
3.6 and 5.8 /im, combined with previous VIA, HST, and Chandra observations. We clearly detect infrared 
emission from the jet, resulting in the most detailed multifrequency data among the jets in lobe-dominated 
quasars. The spectral energy distributions of the jet knots have significant variations along the jet, like the 
archetypal jet in 3C 273. The infrared measurements with IRAC are consistent with the previous idea that 
the jet has two spectral components, namely (1) the low-energy synchrotron spectrum extending from radio to 
infrared, and (2) the high-energy component responsible for the X-ray flux. The optical fluxes may be a mixture 
of the two components. We consider three radiation models for the high-energy component: inverse Compton 
scattering of cosmic microwave background (CMB) photons by radio-emitting electrons in a highly relativistic 
jet, synchrotron radiation by a second distinct electron population, and synchrotron radiation by ultra high 
energy protons. Each hypothesis leads to important insights into and constraints on particle acceleration in the 
jet, as well as the basic physical properties of the jet such as bulk velocity, transporting power, and particle 
contents. 

Subject headings: galaxies: jets — infrared: galaxies — quasars: individual(PKS 1136-135) — radiation 
mechanisms: non-thermal 



1. INTRODUCTION 

Over the last decade, optical and X-ray observations made 
with the Hubble S pace Telescope (e.g., ICrane et ail 119931: 
Bahcall et a l. 1995) and the Chandra X-ray Observatory (e.g., 
Chart aset al.ll2000i) have produced exquisite images of extra- 
galactic kiloparsec-scale jets, completely changing our under- 
standing of their properties. Currently, more than 70 (30) ex- 
tragalactic jets and hotspots are known in the X-ray (optical) 
8 ; all but the few brightest jets were discovered by Chandra 
(HST). 

The origin of the broad-band spectral energy distributions 
(SEDs) of large-scale quasar jets, constructed using HST 
and Chand ra dataware the subject of active debate (fo r re- 
views, see Stawarz 2003J lHarris & Krawczvnskil 120061) . In 
luminous quasars with X-ray jets extending from the quasar 
nucle us out to hundred s of ki lo parsec (e.g.. ISchwartz etal.1 
2000t ISambruna et all |2004t [jorstad & Marscherl 120041 
Marshall et al. 2005; Siemiginowska et al. 2006), the X-ray 
intensity relative to the radio synchrotron flux is generally 
too high to be explained by a synchrotron-self-Compton 
model unless there is a huge deviation from equipartition 
(IChartas et al.l 12000b iKataoka & Stawarzll2005l) . The radio, 
optical, and X-ray fluxes of a jet knot generally trace a 
peaked, inflected broad-band spectrum, which rules out the 
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interpretation of X-rays as due to synchrotron radiation from 
a single population of electrons 9 . An alternative scenario, 
inverse-Compton (IC) scattering of CMB photons by high- 
energy electrons (j e ~ 30) in a highly relativi stic jet with bulk 
Lorentz factor T ~ 10 (the beamed IC model : iTavecchio et all 
l2000tlCelotti. Ghisellini. & Chiabergell2001l) initially seemed 
a more natural way to explain the observed X-ray emis- 
sion, but this process is also not free of problems (e.g., 
lAtoyan & Permed I2004T) . Finally, it is possible that the X- 
rays arise fro m synchrotron r adiation from extremely ener- 
getic protons dAharo nian 2002). Determining which of these 
emission mechanisms produces the observed X-ray jets in 
powerful quasars is a strong motivation for more observations 
of radio-loud quasars, and has resulted in a rapid increase in 
the number of known X-ray and optical jets. 

A new window to explore extragalactic large-scale jets has 
been opened by the Spitzer Space Telescope, which is capa- 
ble of detecting jet infrared emission thanks to t he excellent 



sensit ivity of the Infrared Array Camera (IRAC; Faz io et al 
2004) and Multiband Imaging Photometer (MIPS: IRieke et al 
2004). The first example was the detection of infrared syn- 
chrotron radiation from jet knots in the quasar PKS 0637-752 
with the Spitzer IRAC at wavelengths of 3.6 and 5.8 fim 
(lUchiyama et al. 2005). In terms of the beamed IC model, 
the infrared bandpass is particularly interesting since the bulk- 
Comptonization bump pro duced by cold electrons is ex pected 
to appear in the infrared (Georganopoulos et al. 2005). The 
absence of such features in the PKS 0637-752 jet rules out 
the jet model dynamically do minated by e + e~ pairs in the 
guise of the beamed IC model (Uchiyama et al. 2005). The 
MIPS observations of the j et in Centaurus A (lHardcastle et all 
120061; iBrookes et al.l20 06). and most recently, the results from 
IRAC and MIP S imaging photometry of the jet in M 87 
(IShi et alJ l2006) have been reported, which also demonstrate 

9 It should be noted that, unlike the case of high-power jets, a synchrotron 
interpretation of X-ray em ission has been well est ablished for low-luminosity 
Fanaroff-Riley I jets (e.g., Hardcastle et al. 2001). 
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the powe r of Spitzer to study j et emissions in lower power jets 
(see also lPerlman et al.ll2007l for the M87 jet). 

Now the three Great Observatories collectively offer the 
possibility to identify the radiation mechanisms operating 
in powerful quasar jets. In fact, when c ombined with the 
data f rom the VIA, Hubble, and Chandra dJester et alJl2005l 
2006), the Spitzer observation of the bright quasar 3C 273 
shed new light on the riddle of X-ray jets (Uchiya ma et al.l 
2006). The Spitzer IRAC photometry of the jet knots in 
3C 273 indicated a two-component spectral energy distribu- 
tion: a radio-to-infrared synchrotron component and a sepa- 
rate optical-to-X-ray component. The latter also seems likely 
to be of synchrotron origin, given the similar polarization of 
optical and radio light. The optical polarization, however, has 
not yet been measured with high precision, so this conclusion 
is not yet firm. Perhaps such a double synchrotron scenario is 
applicable to the radiation output from many quasar jets. 

In this paper, we present Spitzer IRAC imaging of the pow- 
erful jet in the luminous quasar PKS 1136 -135. Together 
with data from the VIA, HST, and Chandra (Sambruna et al. 
2002, 2004, 2006a]), our infrared photometry makes the SED 
of the PKS 1136-135 jet the most detailed and best con- 
strained among lobe-dominated quasars. The jet in the 
quasar PKS 1136-135 is reminiscent of the 3C 273 jet 
( Samb runa et al.l l2006al) . demonstrating anti-correlation be- 
tween radio and X-ray brightness, such that the radio inten- 
sity increases toward a hotspot while X-ray flux decreases. 
Applying the beamed IC model to the X-ray emission, this 
has recently been interpreted to imply deceler ation of the jet 
(Sambruna et al. 2006at iTavecchio et al.ll2006l) . Here we ana- 
lyze the multiwavelength jet emission in the light of the dou- 
ble synchrotron scena rio recently outlined for the 3C 273 jet 
(lUchiyamaet al.ll2006l) . The redshift of PKS 1136-135 is 
z = 0.554, so we adopt a luminosity distance of Di = 3.20 Gpc, 
for a concordance cosmology with fi m = 0.27, £!a = 0.73, and 
Hq = 1\ km s" 1 Mpc" 1 . The angular scale of 1" corresponds 
to 6.4 kpc. 

2. OBSERVATIONS AND RESULTS 
2. 1 . Spitzer IRAC observations 

We observed PKS 1136-135 with Spitzer IRAC 
dFazioet al.ll2004l) on 2005 June 10 as part of our Cycle- 1 
General Observer program (Spitzer program ID 3586). We 
used the pair of 3.6 and 5.8 /im arrays, observing the same sky 
simultaneously. The pixel size in both arrays is ~ l."22. The 
point-spread functions (PSFs) are 1 ."66 and lf'88 (FWHM) 
for the 3.6 and 5.8 /im bands, respectively. The photometry 
with I RAC is calibrated to an accuracy of ~ 3% (Rea ch et alJ 
120051) . We obtained a total of 50 frames per IRAC band, 
each with a 30-s frame time. The pipeline process (version 
S 14.0.0) at the Spitzer Science Center yielded 50 calibrated 
images (Basic Calibrated Data). These well-dithered frames 
were combined into a mosaic image with a pixel size of 0."2 
using MOPEX (Ma kovoz & Marleaull2005l) . which removes 
spurious sources such as cosmic rays and moving objects 
based on inter-frame comparisons. Finally, in order to align 
the IRAC image with respect to the VLA image, we set the 
center of the quasar in the IRAC image at the core position in 
the VLA image. In this way IRAC astrometry is estimated to 
be accurate to < Q'.'2. 

The infrared fluxes from the quasar core of PKS 1 136-135 
were measured as 2.2 and 4.0 mJy in the 3.6 and 5.8 /zm 
bands, respectively (Table Q]). These values do not reach the 
saturation limits of the IRAC arrays. Assuming that the in- 
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FIG. 1.— IRAC image of the PKS 1136-135 jet at 3.6 fim in units of 
MJy sr~' with core subtraction. North is up and east is to the left The origin 
of the coordinate is set at the quasar core. A circular region with a radius of 
4" around the core is blanked out, where the core subtraction is incomplete. 
An angular scale 1" corresponds to 6.4 kpc projected size. Superimposed on 
the IRAC image are radio contours on a logarithmic scale from the VLA 8.5 
GHz map. 

frared spectrum has a power law with a = 1 , as consistent with 
the ratio of the 3.6 and 5.8 /im fluxes, the luminosity in the 3— 
10 /im band amounts to Z/s- i Qum ~ 2 .6 x 10 45 ergs s -1 . In 
the model of ISambruna et alJ d2006bl) . the infrared emission 
is assumed to arise from the inner jet, while the optical/UV 
emission originates in the accretion disk. Our infrared mea- 
surement is consistent with that model. If this is the case, the 
infrared fluxes determine the apparent luminosity of the syn- 
chrotron component. 

Inspection of the 3.6/im IRAC images clearly reveals the 
presence of an infrared jet that traces the extended radio jet. 
The IRAC 5.8 /im image also shows the jet emission, al- 
though its quality is worse than the 3.6 /im image. The ex- 
tended PSF of the bright quasar core is significant at the lo- 
cation of the jet, amounting to ~ 10-30% of the jet flux. 
In order to separate the jet infrared emission from the con- 
taminating PSF wings of the core, we subtracted the PSF 
w ings by making use of the PSF templates in a way similar 
to lUchivama et al. (2005). Figure Q] shows the IRAC image 
of PKS 1 136-135 in the 3.6 /im band after subtraction of the 
PSF wings of the quasar core. Infrared counterparts of the jet 
knots are clearly visible in the IRAC image, located ~ 8"— 12" 
from the core (knot B through hotspot HS). A possible coun- 
terpart of knot A is marginally found in the IRAC image, but 
it may be due to our incomplete PSF subtraction, so we re- 
gard the knot A flux only as an upper limit. The hotspot at 
the counter-jet side also suffers from the uncertainties of PSF 
subtraction; we do not find an infrared counterpart and place 
an upper limit of .6 = 4 /iJy at 3.6/im. Also, we do not detect 
extended infrared emissions from the radio lobes. 

To construct the broad-band SEDs of jet knots using mul- 
tifrequency data, we derived infrared flux densities of some 
knot features at 3.6 and 5.8 /im. Given the small separation 
of adjac ent knots, typically ~ 1", we utilized a PSF-fitting 
method (Uchiyama et al. 2006). Specifically, the jet image 
was fitted with a series of IRAC PSFs, fixing the PSF cen- 
ters at the knot positions in the VLA 22 GHz radio image. 
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FIG. 2. — Composite image of the jet in PKS 1 136-135 based on the data 
from Spitzer (shown as red), HST (green dots on knots A and B), and Chan- 
dra (blue). The superposed contours are from the VLA (8.5 GHz), highlight- 
ing the jet radio-emission (see also Fig. [I}. Letters label the main features 
along the jet; knot A closest to the quasar. The multiwavelength image qual- 
itatively illustrates the relative importance of low-energy component (radio- 
infrared) in the outer knots and high-energy one (optical-X-ray) in the inner 
knots. 

However, this procedure did not fully determine knot by knot 
fluxes due mainly to a multiplicity of knots from C to E. 
We then fixed the flux density of knot E at fy.e = 2 /dy and 
/s.8 = 4 /iJy. We derived the systematic errors due to this treat- 
ment by changing the knot E flux in a range of fit = 0-3 /iJy 
and /s.g = 0-6 /Jy. The systematic errors are then added in 
quadrature to statistical ones. The photometric results deter- 
mined in this way are listed in TableQ] where a combined flux 
is reported for knots C, D, and E (referred as knots CDE). 

2.2. Multifrequency Image 

In Figure [2] we present a three-color multifrequency im- 
age based on Spitzer infrared {red), HST optical (green), and 
Chandra X-rays (blue). The VLA 8.5 GHz radio contours 
in a logarithmic scale are overlaid. The infrared photometry 
at 3.6 /im is illustrated as a series of best-fitted PSFs of ev- 
ery knot after artificially shrinking the widths of the PSFs (1" 
FWHM) to restore a resolution similar to the X-rays. At the 
positions of optically bright knots A and B we plot a green dot 
of a ~ Q'. '2 to illustrate clear detection s with the HST STIS 
and ACS dSambruna et alj|2004l2006ah . The Chandra X-ray 
image in 0.3-8 keV band is smoothed with a Gaussian kernel 
of0."4. 

The jet morphology, or brightness pattern along the jet, in 
the infrared appears to be different from emission at shorter 
wavelengths. It is interesting to note that optically-bright 
features (knots A and B) are not particularly prominent in 
the infrared, while an infrared-bright knot D does not emit 
correspondingly strong optical light. The difference in the 
brightness pattern between the infrared and optical suggests 
a dramatic change of the spectral shape in the infrared-to- 
optical band along the jet. This emphasizes the importance 
of measuring infrared fluxes, together with the optical fluxes, 
in quasar jets. 

Putting hotspot HS aside, the jet knots can be divided up 
into two parts: the inner knots (A and B) and the outer knots 
(C, D, and E). The inner knots are bright in both the optical 
and X-rays, and as such they are "high-energy dominated." 
Indeed the inner knots are most luminous in the X-rays, while 
the outer knots seem to have a luminosity peak in the far- 
infrared. To clarify this, we present the SEDs of the jet knots 
in the following section. 

2.3. Spectral Energy Distributions: Two-Component Nature 

We construct the broad-band SEDs of the jet knots using 
the new Spitzer IRAC data ad ded to the VLA, HST, and Chan- 
dra data already presented in Sambruna et al.l d2006a). Once 
again, the Spitzer IRAC provides crucial ph otometric points 
in characterizing the emissions of quasar jets ( lUchivama et alj 



l200ll2006h . 

Table Q] summarizes the flux densities obtained with VLA, 
Spitzer, HST, and Chandra. The radio fluxes measured 
with the VLA at frequencies of 4.9, 8.46, and 22.5 GHz 
are taken from Sambruna et al. (2006a). I nfrar ed photome- 
try with the Spitzer IRAC is described in {32.11 In the opti- 
cal, we list the flux densitie s from the HST ACS in two fil- 
ters dSambruna et al.ll2006ah : F475W (0.48 /im) and F814W 
(0.81 /mi). We employed a some what different m e thod fo r 
X-ray photometry compared with Sambru na et alj ((2006a). 
Specifically, we here define two energy bands for photome- 
try: 1-2 keV and 3-6 keV. Integrating over the entire jet (from 
knot A to hotspot HS) we found 203 and 57 counts in the 1- 
2 and 3-6 bands, respectively (also, 147 counts in the 0.6-1 
keV band). The nominal energy e norm for each band is the 
mean energy weighted by the effective area: e norm = 1 .56 keV 
and 4.35 keV. 

We reana lyzed the Chandra da ta sets that are fully de- 
scribed in Sam bruna et alJ (1200 6a). and constructed flux im- 
ages (in units of photons cm" 2 s -1 pix" 1 ) in the two bands, to 
which aperture photometry was applied. The object apertures 
for knots A/B and hotspot are a circle with a radius of ~ 0"7 
centered on each feature, while that for knots CDE is a box 
enclosing knots C, D, and E. To correct for the interstellar 
absorption of A^h = 3.5 x 10 20 cm" 2 , we multiplied the mea- 
sured fluxes by a factor of 1.04 for the 1-2 keV band. The 
fudge factor slightly depends on the source spectrum such 
that it ranges from 1 .040 to 1 .048 for a power law with pho- 
ton index Y = 1-3. The photometry-based spectra are found 
to be consistent with th e usual spectral analysis presented by 
ISambruna et al.l d2006al) . The photometric method is conve- 
nient when combining lower frequency data. 

Care should be taken in interpreting the X-ray fluxes. First, 
the X-ray emission downstream of knot B does not show well 
discernible knots, and therefore the X-ray flux of knots CDE 
may not be related directly to the lower frequency fluxes. Sec- 
ond, in knot A ther e is an offset of ~ 0."4 b etween the soft and 
hard X-ray peaks (Samb runa et al.ll2006al) . The X-ray spec- 
trum derived for knot A may be contaminated by some unre- 
lated emission, which makes the spectral shape inconsistent 
with a power law. 

In Fig. [3] we present the SEDs from radio to X-rays for 
knots A, B, CDE, and hotspot HS. We note that the infrared 
fluxes we have measured with the Spitzer IRAC fill central 
points in the radio-to-X-ray SEDs of the jet knots, setting an 
important constraint on models of the broad-band emission. 
The inner knots (A and B) radiate most strongly in the X- 
ray. Also, significant (and probably flatter) optical emission 
is observed exclusively from the inner knots. The possible dif- 
ference in slope s of the optical continu um between the inner 
and outer knots (Sambruna et al. 2006a) suggest different ori- 
gins; the optical emission in the inner knots may belong to the 
same spectral co mponent responsible for the X-ray emission 
(as suggested by Sambr una et alJl2006aL for knot A), while 
in the outer knots the optical emission is related to the radio 
synchrotron component. 

These spectr al characteristics are an alogous to those of the 
jet in 3C 273 (lUchivama et al.ll2006l) . for which we identi- 
fied two spectral components: (1) the low-energy synchrotron 
spectrum extending from radio to infrared, and (2) the high- 
energy component arising in the optical and smoothly con- 
necting to the X-ray flux. We argued that the second compo- 
nent is likely to be of synchrotron origin as well, thus form- 
ing double-synchrotron spectra, because of the similarity of 
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TABLE 1 

Flux Densities of the PKS 1 136-135 Core and Its Jet Knots 



Flux Density, f v 



Region /s GHz k,i GHz /22 GHz As Aim fi.b fim /o.81 fjm /o.48 pm /l.56kcV ./*4.35 keV 

(mjy) (mJy) (mJy) (/xjy) (p]y) (^tjy) (p]y) (nJy) (njy) 



Core 

A 

B 
CDE 

HS 



463 

11 ±1.6 
192 ±19 
200 ± 20 



3.8 ±0.4 
9.3 ±0.9 
116±12 
1 19 ± 12 



3.2±0.6 
47.5 ±9.5 
43.9±8.8 



4.0 x 10 J 
< 5 
4.8 ±2.2 
20±4 
7.4 ±2.7 



2.2 x 10 3 
<4 
3.8± 1.9 

9.5 ±2.0 

4.6 ±2.0 



0.23 ±0.07 
0.33 ±0.09 
< 0.46 
0.24 ±0.05 



1.5 x 10 3 
0.17 ±0.03 
0.22 ±0.03 
<0.23 
< 0.08 



170 

0.77 ± 0.22 
2.0 ±0.34 
1.2 ±0.26 
<0.16 



0.72 ±0.29 
0.72 ±0.24 
0.67 ±0.25 



NOTE. — The radio and optical data (extinction-corrected) for jet knots are taken from Sambruna et al. 1 2006a). The errors of 
optical fluxes are translated into those at a 90% level. The infrared fluxes of the jet knots are determined through PSF-fitting and 
quoted errors are dominated by their systematic errors. The X-ray fluxes listed with statistical errors (at a 90% confident level) are 
corrected for Galact ic absorption using Nu = 3.5 X 10 20 cnT 2 . Except for the infrared, the data of the core fluxes are taken from 
Gambill et al. 1 2003) assuming a opt = 0.5 and a\ = 0.7. 



TABLE 2 

Phenomenolooical Model Parameters of PKS 1136-135 Jet 



Parameter 


knot A 


knot B 


knots CDE 


hotspot HS 


Ctl 


0.65 


0.65 


0.85 


1.0 


y\ (Hz) 


< 3 X 10 13 


2x 10 13 


3 X 10 13 


5 X 10 13 


Li (ergs -1 ) 


< 1.3 X 10 43 


3.1 x 10 43 


1.7 X 10 44 


1.2 X 10 44 


Oil 


0.75 


0.70 


0.50 




U (erg s" 1 ) 


5.0 x 10 43 


8.0 x 10 43 


7.3 X 10 43 





NOTE. — Equation {TJ is adopted to characterize the SEDs. 

polarization between optical and radio emission (though the 
degrees of optical polarization are not established yet). 

Given the similarity to 3C 273, we model the radio-to-X-ray 
SEDs phenomenologically by the following function, namely, 
a double power-law with an exponential cutoff: 



fu = ^2 K ' v °" exp 



1/2' 



(1) 



The first term (i = 1) of the right-hand-side accounts for the 
low-energy spectrum and the second term (i = 2) describes 
the high-energy part. As there are no firm indication of 
spectral steepening at X-rays, we set an arbitrary cutoff at 
i>2 = 1 x 10 Hz for the second component. For the inverse- 
Compton model described later ( 93.11 see Fig. |5), however, 
the high-energy cutoff is located far beyond the X-ray domain, 
reaching very high-energy gamma-rays. 

Table [2] lists the spectral parameters that well describe the 
knot SEDs shown in Fig. [3] Instead of the normalization, k,, 
we present the total apparent luminosity (without taking ac- 
count of relativistic beaming) defined as L = AnD\§ f v dv for 
each component, where the low-energy boundary of the in- 
tegral is set at 10 9 Hz and 10 13 Hz for the first and second 
component, respectively. Knot B, the most X-ray luminous 
one, has a luminosity of L2 ~ 8 x 10 43 erg s" 1 for the high- 
energy component, which indeed exceeds the luminosity of 
the low-energy component. 

The spectral index of the low-energy component is in the 
range of ot\ ~ 0.65-1.0, which is determined by the radio 
spectra, as listed in Table [2] The cutoff frequency of the 
low-energy component, as constrained by the IRAC fluxes, 
is v\ = (2-5) x 10 13 Hz in the cases of knots B, CDE, and 
hotspot HS; v\ < 3 x 10 13 Hz for knot A. The second spectral 
index is found to be «2 — 0.5-0.75. 

3. INTERPRETING THE HIGH-ENERGY COMPONENT 



While the low-energy emission from quasar jets is un- 
doubtedly of synchrotron origin, the production mech- 
anism (s) of strong X-rays rema ins a matter of debate 
(e.g., lHarris & Krawc zynski 2006). Three radiation models 
have been proposed to explain the high-energy component: 
inverse-Compton scattering by radio-emitting electrons, syn- 
chrotron radiation by a second electron population, and syn- 
chrotron radiation by energetic protons . In this section we 
consider each of these in turn. 

There are reasons to think t hat jets are relativistic on large 
scales (Garrington et al. 1988) and thus beaming may be im- 
portant. The Doppler beaming factor is defined as 5 = [T{1 - 
ficosO)]' 1 where f3c is the velocity of the jet, T = (1 - /3 2 )" 1/2 
is the bulk Lorentz factor of the jet, and 9 is the observing 
angle with respect to the jet direction. A one-sided jet with 
two-sided hotspots in this source is a sign of substantial rel- 
ativistic beaming. The Doppler factor is, however, not well 
constrained for this jet. Throughout this section, all phys- 
ical quantities are referred to a jet co-moving frame unless 
otherwise specified. The exceptions to this are the direction 
and velocity of the jet itself. It is assumed that a relativisti- 
cally moving blob (as noted by knot or hotspot) has spherical 
geometry with a radius of r = 1 kpc (in the jet frame of ref- 
erence) 11 occupied homogeneously with relativistic particles 
and magnetic fields ("one-zone" model with a filling factor of 
unity). 

Before we discuss the high-energy emission component of 
the jet knots (like knot B), we derive a lower limit on the 
magnetic field strength in hotspot HS, based on the absence 
of the high-energy emission. We do not introduce beam- 
ing effects for hotspot HS (T = 6 = 1). Synchrotron pho- 
tons are Compton upscattered by synchrotron-emitting elec- 
trons themselves, namely "synchrotron self-Compton" (SSC), 
which is known to be a primary component for the high- 
en ergy emission in bright hotspo ts such as those in Cygnus 
A (lHarris. Cari lli. & Perley 1994). In fact, since the average 
energ y density of synchrotron photons (see lBand & Grindlavl 
119851) in hotspot HS, U sya ~ 9L syn /(167rr 2 c) ~ 100 eV cm" 3 , 
largely exceeds the CMB energy density of C/cmb = 0.26(1 + 
z) 4 eV cm" 3 ~ 1.5 eV cm" 3 , the SSC flux well exceeds the 
IC/CMB flux in the hotspot of PKS 1136-135. We found 



10 Yet another possibility includes the model of Dermer & Atovan (2002), 
which invokes spectral hardening due to inverse-Compton cooling in the 
Klein-Nishina regime. 

11 The size of the emission region is not certain, as the knots are unre- 
solved. The knot size is assumed to be similar to the 3C 273 jet. 
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: (a) knot A 




FIG. 3.— Broadband radio-to-X-ray SEDs for knots A, B, CDE, and 
hostspot HS in the powerful quasar 1 136—135. IRAC 3.6 and 5.8 fim fluxes 
are shown as open circles. The data are su mmarized in TableQ] In addition, a 
flux density of 12.1 ± 1.8 mJy at 1.5 GHz iSambruna et al. 2006a) is plotted 
in the SED of knot A. The two-component phenomenological model given as 
equation (TJ is used to describe the data. 

a lower limit on magnetic field as B > 0.2 mG in order not 
to violate the X-ray upper limit by the SSC emission. With 
an equipartition magnetic filed of B eq = 0.66 mG where the 
energy density of relativistic electrons is equal to that of mag- 
netic fields, the one-zone SSC flux is an order of magnitude 
lower than the X-ray upper limit. 

3.1. Inverse-Compton Scattering by Ultra-Relativistic Jet 

We first discuss the second, X-ray-dominated component 
of the SED seen in the jet knots (knot B, in particular) 
in the framewor k of the b eamed inverse-Compton model 
dTavecchio et alj l2000t ICelotti. Ghisellini. & Chiabergd 



1200 lb . which has been advocated in the previous pape rs for 
this jet (ISambruna et al.l l2006ah ITavecchio etail 120061) . In 
this model, it is assumed that the jet has a highly relativistic 
bulk velocity, with a Lorentz factor T ~ 10 out to distances 
of hundreds of kiloparsecs, so that the CMB field seen by 
electrons is enhanced by a factor of T 2 in the jet comoving 
frame. The amplified CMB can be Compton up-scattered by 
high-energy electrons of 7 e ~ 30 into the observed X-rays. 
This model automatically requires that the jet direction is 
close to the line-of-sight, 9 ~ T" 1 . (Such a condition, and 
therefore the beamed IC model itself, would not be favorable 
for a lobe-dominated quasar like PKS 1136-135. Later in 
this section, we shall briefly discuss this issue.) 

It is interesting to test the beamed IC/CMB hypothesis, 
particularly for the PKS 1136-135 jet. The jet exhibits the 
evolution of multiwavelength emission along the jet, namely 
the increase of radio brightness accompanied by the decrease 
of X-ray brightness towards dow nstream. As detailed in 
Georganopoulo s~& Kazanasl (120041) . such a behavior can be 
interpreted, within the framework of t he beamed IC radiation , 
as being due to decele ration of the jet (Sambruna et al. 2006a; 
ITavecchio et al J 120061) . Specifically. ITavecchio et al.l (2006) 
concluded that the bulk Lorentz factor of the jet decreases 
from r ~ 7 (knot B) to T ~ 2 (knot E). If the beamed IC sce- 
nario is confirmed, we gain a new probe of the flow structure 
of the large-scale jets. 

Let us apply a synchro tron plus beamed IC/CMB model 
of lTavecchioeT al. (2000) to the broadband SED of knot B. 
We also include SSC calculation for completeness. The both 
models adopt "one-zone" emission, assuming an emitting 
blob with homogeneously filled with relativistic particles and 
magnetic fields. It is interesting to note that in the case of 
no-beaming (r = 6 = 1), again, the average energy density of 
synchrotron photons in the knot, £/ syn ~ 12 eV cm -3 , exceeds 
the CMB energy density. The energy distribution of radiating 
electrons is assumed to be of the form N(sy) = kj~ s , having a 
low-energy cutoff at 7 m ; n = 20 and a high-energy exponential 
cutoff at 7 max , where 7 denotes the electron's Lorentz factor. 
The number index of electrons is set to be s = 2a + 1 =2.4 
based on a typical radio index in this jet (a ~ 0.7), and the 
high-energy cutoff is determined to be 7 max = 4 x 10 5 to pro- 
duce the infrared flux. The magnetic field strength is adopted 
as B = B e q,5=i^ _1 5 where B eq .,5 = i = 175 /iG; this relation en- 
sures equipartition between relativistic electrons and mag- 
netic fields, and keeps the peak frequency of the synchrotron 
radiation, j/ max oc SBj^^. The bulk Lorentz factor itself does 
not explicitly affect the discussion made here and it is adopted 
to be T= 15. 

Figure g] shows the IC/CMB and SSC X-ray fluxes as a 
function of 5, given the radio flux observed for knot B. Since 
the normalization of the electron distribution scales as k oc S~ 2 
and the beaming pattern of the IC/CMB emission follows <5 3+ J 
dDermerlll995HGeorganopoulos. Kirk. & Mastichi adis 2001), 
the IC/CMB flux scales as kS 3+s oc S l+S = <5 3 4 . The SSC flux 
scales as k 2 B 1+a S 3+a oc <T 2 , where a = (s- 1)/2 = 0.7. As 
such, if the jet is heavily beamed, 6 ^> 1, the IC/CMB flux 
is dominant. It reaches the observed X-ray flux at 5 = 18.6; 
the corresponding SED is shown in the top panel of Fig. [5] 
The SSC becomes a main component for S < 1 .4. Even if we 
invoke a particle-dominated case of U e /Us = 5000, the SSC 
model requires a significant de-beaming of 6 = 0.13 to explain 
the observed X-rays, which is quite unlikely as the on-axis 
version of such jet emission becomes unacceptably luminous. 

The radio-emitting electrons emit GeV 7-rays through the 
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FIG. 4.— Predicted flux of IC/CMB and SSC at 1 keV as a function of the 
beaming factor (<5) based on the observed radio flux of knot B. The observed 
X-ray flux is 7 X 10 s JyHz (a horizontal dotted line). The IC/CMB emission 
is dominant for a large <5, while SSC dominates the X-ray emission for de- 
beaming cases (<5 < 1). A dashed line shows the SSC prediction with the 
assumption of a significant deviation from equipartition (U e /Ub = 5000). 

IC/CMB process (Pig. [5j. The total radiative output is dom- 
inated by the multi-GeV 7-rays and their predicted flux just 
reaches the sensitivity offered by the upcoming GLAST satel- 
lite (see lGeorganopoulos et al.ll2006l for the case of 3C 273). 
The quasa r core is expected to sh ow a similar level of 7-ray 
emission ( Sam bruiia et alj|2006bl) : steadiness and hard spec- 
trum are signatures of the large-scale emission. The observa- 
tion of PKS 1 136-135 (or quasar jets in general) with GLAST 
may be able to test the IC/CMB hypothesis. Note that the sub- 
TeV domains are subject to intergalactic absorption, which 
makes it difficult to investigate the jet emission with future 
ground-based Cherenkov telescopes. 

The Doppler factor of S ~ 19 required by the IC/CMB 
model implies an uncomfortably smal l angle betw e en the 
jet and the line-of-sight, 9 < 4°. iTavecchio et alJ (2006) 
also derived the maximum angle permitted as 9 max = 3.°8 for 
this jet. The small angle would make the jet quite long, 
> (77 / sin^max) kpc =1.1 Mpc, and the total source extent 
would be as large as the largest quasar. The core-to-lobe 
flux ratio at 5 GHz, R, is an indicator of jet orientation. In 
the case of 1136-135, the cor e-to-lobe ratio of R = 0.29 is 
obtained (Samb ruiia et al .1120041) . t hus designated as a lobe- 
domin ated quasar. According to lLaurent-M uehleise n et alJ 
d!997l) . a sample of radio-loud quasars with a median value 
of R = 0.52 have an average angle to the line-of-sight of 
(9) ~ 30° with simple beaming models. Also, the mean an- 
le for steep-spectrum r adio quasars is estimated as (9) ~ 30° 
Urry & Padovanil 19951) . In this respect, a very small viewing 
angle 9 < 4° for a lobe-dominated quasar like PKS 1 136-135 
is not favored in the light of simple un ification schemes 
of rad io-loud quasars and radio galaxie s dUrrv & Padovan 
9951): a similar statement was made by Kataoka & Stawar: 
(2005). However, it should be kept in mind that Chandra 



detected jets can be biased toward smaller jet angles, which 
alleviates a problem of the jet angle. We need a more direct 
means to distinguish currently proposed models (see 33.41 ). 

3.2. Synchrotron Radiation by Second Electron Population 

The large-scale jet in PKS 1136-135 is in many re- 
spects similar to the we ll-known jet of quasar 3C 273. As 
Sambruna et al.l d2006al) have pointed out, in both jets, the 
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FIG. 5. — (Top) Synchrotron+IC/CMB modeling (solid and dashed curves) 
for knot B. The SED points are same as Fig. [5] A large beaming factor of 
6 = 18.6 is required to explain the X-ray flux. For illustrative purpose, it is 
also shown that if instead we invoke de-beaming (<5 = 0.1: a dotted curve), 
the SSC luminosity becomes comparable to the synchrotron one. The ex- 
pected sensitivity of the GLAST satellite (1 year 5cr) is marked. (Bottom) 
Proton synchrotron model for the high-energy emission with moderate beam- 
ing r = 8 = 3 (6 ~ 19°). The low-energy emission is ascribed to electron 
synchrotron radiation. The dotted curve shows the electron component alone, 
while the solid line presents synchrotron radiation produced by both protons 
and electrons. Assuming equipartition between protons and magnetic fields, 
the magnetic field strength is B c q = 1.7 mG. The production rate of electrons 
in terms of the energy content is about 5% of that of protons. 

radio emission brightens monotonically towards the terminal 
hotspot, while the upstream knots are brighte r in X-rays. The 
two-c omponent SEDs are also similar (cf, lUchiyama et al.1 
2006). In the case of 3C 273, since both the radio and 
optical emission are linearly polarized to a similar degree 
and in the same direction, it seems reasonable that the high- 
energy component, which contributes more than half the op- 
tical flux, should also arise from the synchrotron process 
dUchiv ama et al. 2006). Therefore we discuss the optical-to- 
X-ray component of 1 136-135 in terms of synchrotron radi- 
ation produced by a second population of high-energy elec- 
trons. The suggestion of the presence of a second synchrotron 
component has been made in previous studies of other large- 
scale jets (e.g.. iJorstad & Marschei1l2004t [Harris et al.ll2004t 
ICheung et al.ll2006l) . 

Let the energy distribution of accelerated electrons be 
Q(E e )dE e oc E~ Sl dE e for an energy interval of interest. The 
power-law slope of electrons is derived from the radiation 
slope as 52 = 2a?2 [see Eq. (Q]l] in a synchrotron-cooling 
regime as appropriate for optical and X-ray emission; the 
cooled electrons responsible for the observed radiation have 
a stee per number index of S2 + 1 (see e.g., iKataoka et al.l 
120031) . Then the spectral index of the high-energy compo- 
nent, «2 — 0.5-0.75 corresponds to S2 — 1.0-1.5. However, 
within the framework of diffusive shock acceleration, the 
hardest possible electron distribution of acceleration would be 



Infrared Study of PKS 1 136-135 Jet 



7 



Q(E e ) oc E~ 15 , in other words, S2 > 1.5 (see Mal kov & Drurvl 

2001) . The inferred index 52 — 1.0-1.5 would violate such a 
theoretical limit. This brings into question the idea that the 
second electron distribution is formed through the diffusive 
shock acceleration or that there is a second electron distribu- 
tion at all. Note that the jet in 3C 273 has spectral index of 
«2 — 0.7-0.8 corresponding to 52 — 1.4-1.6, which was con- 
sidered to be compa tible with the shock acceleration theory 
dUchivama et al.ll2006l) . 

The second synchrotron component in the PKS 1136-135 
jet may instead be due to turbulent acceleration operating 
in the shear layers (Ostrowski 2000; Staw arz & Ostrowskl 

2002) . The coexistence of two distinct types of accelera- 
tion, shock and turbulent acceleration, in the knots would 
naturally give rise to a double-synchrotron spectrum. More 
importantly, unlike the shock acceleration, a hard spectrum, 
52 < 1-5, can be expected to form in the case of turbulent 
acceleration (i.e., second-order Fermi acceleration). In fully 
turbulent shear layers, the time scale of turbulent accelera- 
tion can be estimated as f acc ~ r g cv~^, where r g = E e /{eB) is 
an electron gyroradius and va denotes the Alfven velocity. 
The turbulent acceleration has to compete against synchrotron 
losses, which is presumed to dominate over IC losses, with a 
timescale of 



(4/3W/B7f 



(2) 



where aj is the Thomson cross section, and Ub = B 2 /(8tt) is 
the energy density of the magnetic field. By equating t acc = 
fsyn — namely, balancing the acceleration and synchrotron loss 
rates — one obtains the maximum attainable energy limited by 
synchrotron losses: 



6x 10 



15 V A 



1 



eV, 



(3) 



where B-4 = B/(1Q 4 G). To account for the observed 
X-rays with synchrotron radiation, one needs E ermx > 
30 (B- 4 ST - 5 TeV fsee IUchivama1it~aT1 [2006). thus requiring 
va > 5 x 10~ 3 c. This condition is reasonable with typical jet 
parameters. 

3.3. Proton Synchrotron Radiation 

Finally, we argue the possibility that the optical-to-X-ray 
emission may be due to synchr otron radiation by very high en- 
ergy protons (I Aharonianl2 002 ) . This model requires that pro- 
tons in the jet are somehow accelerated to very high energies, 
E p ~ 10 18 eV or more. Also, the magnetic field strength must 
be of the order B ~ mG, so that energy equipartition between 
the relativistic protons and magnetic fields can be (roughly) 
realized. The knots and hotspots in the relativistic jets of pow- 
erful quasars and radio galaxies are indeed one of a few po- 
tential sites o f cosmic-ray acceleration up to E p ~ 10 20 eV 
dHiilaslll984l) . For example, the turbulent acceleration in the 
shear la yer may be able to accelerate such ultra high energy 
protons dOirowski 2003). 

The characteristic frequency and cooling time of proton 
synchrotron radiation can be written as 



and 



^p-syn - 



'p-syn 



2.1 x 10 1 



1.4 x 10 s 



B 

mG 



B 

mG 



10 18 eV 



10 18 eV 



Hz, 



yr. 



(4) 



(5) 



respectively (lAharonianl 12002). (Note that the synchrotron 
cooling time is much shorter than the photo-meson cooling 
time at energies relevant here.) A total energy content of the 
magnetic field B = 1 mG integrated over a spherical knot with 
a radius of r = 1 kpc amounts to Wb = 4.9 x 10 57 erg. If we 
assume equipartition between radiating protons and magnetic 
fields, W p = Wb, the X-ray luminosity of synchrotron radia- 
tion by ultra-high-energy protons of E p ~ 10 18 eV can be es- 
timated roughly as 



( p-syn 



4.9 x 10 57 erg 
1.4 x 10 8 yr 



V 



1.1 x 10 42 77 ergs" 1 , (6) 



where 77 ~ 0.1 denotes the fraction of kinetic energy of pro- 
tons responsible for the X-ray emission. Relativistic beaming 
is not taken into account in this estimate. Roughly speaking, 
the luminosity scales with B and 5 such that L x oc (B5) 4 . The 
expected luminosity agrees with typical (apparent) luminosity 
of the jet knots of quasars, L x ~ 10 43 erg s -1 , if BS ^ 6 mG. 
Thus, very high energy protons rather than a second popula- 
tion of electrons, may explain the second synchrotron compo- 
nent extending from the optical to X-rays provided that both 
ultra-high-energy protons of E p ~ 10 18 eV and a magnetic 
field of B ~ mG are present in the jet knots. Under the proton- 
synchrotron hypothesis, the low-energy synchrotron emission 
may be accounted for by accompanying electrons. 

To be specific, we modeled the high-energy spectral com- 
ponent of knot B with proton synchrotron radiation. We adopt 
mild beaming of S = T = 3 (6 ~ 20°). The low-energy spec- 
tral component was simultaneously modeled by electron syn- 
chrotron radiation with taking account of the effects of signif- 
icant synchrotron cooling. It is assumed that during Af = 10 7 
years, protons and electrons are injected continuously into the 
emission volume (r = 1 kpc) with the distributions character- 
ized by number index s p and s e and by the maximum energy 
E p . max and E e , max : Q(E) = kE~ s exp(-£ /E max ). In the bottom 
panel of Fig. [5] we present the broadband SED reproduced by 
the proton synchrotron model with the equipartition magnetic 
field of B = 1.7 mG, the index of the power-law distribution of 
accelerated particles s p = s e = 1.8, and the maximum energy 
of is^max = 4.7 x 10 17 eV and £ c>m ax = 1 x 10 11 eV. The injec- 
tion rate of electrons, in terms of energy content J Q(E)dE, 
is about 5% of that of protons. The injected electrons suffer 
from severe synchrotron cooling. Even radio-emitting elec- 
trons have a cooling time of only ~ 10 4 years, which requires 
in-situ acceleration of electrons. The infrared-emitting elec- 
trons have a lifetime of ~ 10 2 years, implying that each knot 
is a currently (within 10 2 years) active site of particle acceler- 
ation. 

The energy density of protons U p , which is equal to 
U B , can be calculated by U p ~ (At/V) J Q p (E)dE ~ 7.4 x 
10 4 eV cm" 3 with V = 4irr 3 /3, since synchrotron cooling of 
protons is not effective. The kinetic power of the jet, es- 
timated as Lj et ~ irr 2 T 2 {3c(U p + Ub), amounts to Lj et ~ 2 x 
10 48 erg s" 1 . This estimate is not sensitive to the choice of 
r as long as the equipartition condition is fulfilled. A large 
power has to be carried by the protons and magnetic fields in 
the jet to produce the X-ray emission via a proton-synchrotron 

lat power s 
10 120021) 

relation to the width of the H/3 line and the optical continuum 
luminosity. The black hole mass corresponds to the Edding- 
ton luminosity of Lndd — 6 x 10 46 erg s . Therefore, the pro- 
ton synchrotron model requires the jet to be super-Eddington, 



process. The black h ole mass that power s the jet is estimated 
to be 4.6 x 10 8 M Q (Oshl ack etail. 12002) using the empirical 
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Ljet ^> ^Edd- This issue would pose a problem for the pro- 
ton synchrotron model, though it may be alleviated by assum- 
ing that the knots present the locations of power peaks due to 
modulated activity of the central engine. 

3.4. Possible Discriminator of the Radiation Models: 
Polarization of Optical Light 

We emphasize here that optical polarimetry can be an effec- 
tive way of discriminating the radiation models responsible 
for the optical-to-X-ray emission of the jet in PKS 1 136-135, 
and of quasar jets in general. In the beamed IC interpretation, 
if the optical fluxes belong to the IC component, the optical 
and X-ray emission are due to Compton up-scattering off the 
amplified CMB by high-energy electrons of r y e ~ 3 (optical) 
and 7e ~ 100 (X-ray). Unlike in the case of synchrotron mod- 
els, the X-rays are expected to be unpolarized and the optical 
light is nearly unpolarized at most a few percent of polariza- 
tion (Uchiyama & Coppi, in preparation) 12 . Precise polar- 
ization measurements in the optical can in principle verify (or 
discard) the beamed IC model. Unfortunately, there have been 
no useful polarization observations of quasar jets with HST so 
far. Only for 3C 273, early HST polarimetry of the jet was 
done but with the pre -COSTAR FOC and low significances 
dThomson et alJfl993|). As was recentl y performed for nearby 
radio galaxies by lPerlman et alJ(l2006l) . new, deep polarimetry 
of quasar jets on large-scales is highly encouraged, although 
the opt ical emission from the jet knots is generally faint, say 
< /Jy (ISambruna et alJl2004h . 

4. SUMMARY 

We have detected the infrared emission from the kiloparsec- 
scale jet in the quasar PKS 1 136-135 with the Spitzer Space 
Telescope. Using the new Spitzer data together with deep 
Chan dra, Hubble, and multi -frequency VLA observations 
from ISambruna et al.l ((2006a), we construct the broadband 
spectral energy distributions along the large-scale jet. The 
SEDs of the jet knots are comprised of two components: 
radio-to-infrared synchrotron emission and a separate high- 
energy component responsible for X-rays, which may extend 

12 Compton up-scattering by a highly relativistic jet in which a popula- 
tion of cold electrons with y e ~ 1 (in the jet frame) exists, namely "bulk 
Comptonization", can yield in principle a high degree of polarization for a jet 



down to the optical. In total, each component has a similar 
apparent luminosity of the order of L ra( j ~ 10 44 erg s" 1 . 

The origin of the high-energy component has significant 
implications to the properties of particle acceleration in a rela- 
tivistic jet as well as jet dynamics. We consider three radiation 
models that have been proposed to explain the high-energy 
component: inverse-Compton scattering in a highly relativis- 
tic jet, synchrotron radiation by a second electron population, 
and synchrotron radiation by highly energetic protons. 

In terms of the double synchrotron scenario, a flat spec- 
trum of the second synchrotron component, cti_ ~ 0.5, would 
turn down diffusive shock acceleration as the mechanism of 
electron acceleration to very high energies, > 30 TeV. In- 
stead turbulent acceleration in the shear layers may explain 
the second population of electrons giving rise to the flat spec- 
trum. In the beamed IC interpretation, a jet orientation with 
respect to the observer of 9 < 4° is required, which is not 
quite favorable given the small core-to-lobe ratio of this ob- 
ject. Future optical polarimetry or 7-ray observations with 
GLAST may be able to test the beamed IC model. The pro- 
ton synchrotron model requires a relatively large magnetic 
field of ~ 1.7 mG, and consequently large kinetic power of 
Ljet — 2 x 10 48 erg s" 1 ~ 30LEdd, with efficient acceleration of 
protons to E p max ~ 5 x 10 17 eV. If this is the case, one can 
probe particle acceleration at very high energies, higher than 
any other known objects in the Universe. 
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ana le 9 ~ r-' (Begel man & Sikora 1987). For 7^ of a few as relevant here, 
however, the polarization is largely suppressed. 
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